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I.  INTRODUCTION 


In  recent  years,  a number  of  methods  have  been  developed  In  order  to 
compute  laminar  and  turbulent  three-dimensional  boundary  layers.  Although 
some  are  developed  in  a general  way,  most  are  developed  to  treat  either  of 
two  specific  classes  of  problems,  i.e.,  flows  on  wings,  or  flows  past  bodies 
or  revolution,  Cebeci  et  al.  [1]  developed  a general  method  for  calculating 
the  three-dimensional,  compressible  boundary  layer  on  an  arbitrary  wing  using 
the  eddy-viscosity  concept.  Similar  calculations  have  recently  been  presented 
by  McLean  [2]  and  by  Kordulla  [3].  Bradshaw  et  al.  [4]  simplified  the  three- 
dimensional  equations  on  a swept  tapered  wing  by  assuming  conical  similarity, 
and  used  a higher  order  turbulence  model  to  solve  for  the  flow.  The  history 
of  flows  past  bodies  of  revolution  is  more  extensive,  dating  back  to  the  first 
studies  by  Moore  [5]  dealing  with  flow  past  a cone.  Further  studies  by  Cooke 
[6],  Boericke  [7]  and  others  integrated  the  conically  similar  quasi-two- 
dimensional  equations  around  the  cone.  The  spinning  cone  and  its  attendant 
problems  were  investigated  by  Dwyer  [8]  and  Kitchens  et  al.  [9].  The  diffi- 
culties encountered  in  all  the  previous  works  due  to  behavior  near  the  lee 
side  of  the  cone,  or  crossflow  reversal,  were  effectively  overcome  by  Lin  and 
Rubin  [10],  [11]  who  developed  an  augmented  set  of  boundary-layer  equations 
for  treating  problems  of  this  kind. 

The  flow  past  a more  general  body  of  revolution,  a prolate  spheroid,  has 
been  extensively  studied  by  Wang  [12]  - [14].  Although  the  cone  is  a special 
case  of  a general  body  of  revolution,  for  the  prolate  spheroid  at  least,  some 
of  the  more  troublesome  aspects  of  the  cone  flow  do  not  appear.  The  nonexist- 
ence of  solutions  on  the  lee  plane,  found  by  Moore  [5]  and  further  investigated 
by  Roux  [15]  and  Murdock  [16],  does  not  seem  to  apply  to  the  prolate  spheroid. 
The  crossflow  separation,  which  causes  the  quasi -two-dimensional  conical  flow 
to  experience  a square-root  singularity,  does  not  seem  to  exert  such  an 
effect  on  the  equations  governing  the  prolate-spheroid  flow.  This  has  allowed 
Wang  to  carefully  map  the  attached  flow  on  the  prolate  spheroid  at  incidence 
by  integrating  from  either  the  windward  or  leeward  synmetry  plane,  something 
which  cannot  generally  be  done  on  a cone  using  standard  boundary-layer 
procedures . 
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The  general  characteristics  of  the  viscous  flow  past  a body  of  revolu- 
tion can  be  broadly  described  as  follows.  Depending  on  the  incidence  of  the 
body,  the  flow  can  be  thought  to  consist  of  three  overlapping,  and  not-at-all 
distinct,  regimes:  (a)  low  incidence,  where  no  separation  occurs  and  three- 
dimensional  boundary-layer  theory  applies;  (b)  low  to  moderate  incidence, 
where  the  cross  flow  around  the  body  separates  and  may  possibly  produce  stream- 
wise  vortices  which,  however,  still  remain  within  the  viscous  layer;  (c)  higher 
incidences,  where  vortices  detach  from  the  body  and  alter  the  prevailing 
inviscid  flow.  Case  (c)  requires  the  complete  (and  pf>ssibly  unsteady)  Navier- 
Stokes  equations,  but  case  (b)  can  be  described  by  a simpler  set  of  equations 
called  here  the  Parabolic-Elliptic  Boundary-Layer  (PEBL)  equations. 

Under  the  research  sponsored  by  contract  N60921-76-C-0089,  aspects  of 
this  problem  were  investigated  at  Douglas  during  the  past  year,  with  the 
emphasis  being  towards  the  development  of  a general  computer  program  for  cal- 
culating the  boundary  layer  on  an  arbitrary  body  of  revolution  using  the  clas- 
sical three-dimensional  boundary-layer  equations,  i.e.,  case  (a).  A parallel 
effort  used  the  PEBL  approach  in  order  to  extend  the  calculations  to  case  (b). 
This  report  presents  the  results  of  the  first  phase  of  the  above  contract. 
Calculations  were  performed  on  a simple  body  of  revolution,  the  prolate  spher- 
oid. For  the  case  of  the  standard  three-dimensional  boundary-layer  equations, 
the  prolate  spheroid  produces  a region  of  reverse  crossflow,  due  to  an  adverse 
pressure  gradient,  but  without  a change  of  sign  of  the  cross-flow  edge  velocity. 
Thus,  the  method  used  in  [1],  where  the  cross-flow  integration  direction  is 
changed  when  the  edge  velocity  changes  sign,  cannot  be  used,  and  the  behavior 
of  the  finite-difference  solution  technique  becomes  apparent.  The  PEBL 
approach  is  used  for  the  first  time  on  a general  body  with  variable  metric 
coefficients  (not  simply  a cone)  and  the  sensitivity  of  the  procedure  to 
streamwise  separation  can  be  assessed. 

A supplementary  effort  dealt  with  the  incorporation  of  a new  scaling  law 
into  the  quasi -two-dimensional  conical  boundary-layer  equations  for  the 
turbulent  flow  on  a cone  [17].  A suggestion  made  by  Bradshaw  [4],  concerning 
the  choice  of  similarity  variable  to  be  used  in  the  computation,  was  followed, 
and  produced  noticeably  better  comparisons  with  the  experimental  data.  These, 
and  other  results  from  the  various  programs  will  be  described  in  the  succeeding 
sections. 
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II.  GOVERNING  EQUATIONS 


2.1  The  Boundary-Layer  Equations  for  a General  Body  of  Revolution 

In  this  report  we  consider  the  solution  of  more  general  three-dimensional 
boundary-layer  equations  than  those  discussed  in  our  previous  studies  [1],[18], 
[19].  The  generalization  is  due  to  the  inclusion  of  crossflow  diffusion  terms.  A 
detailed  derivation  is  given  in  Appendix  A.  These  generalized  equations  allow 
for  the  computation  of  negative  crossflow  due  to  the  elliptic  nature  of  the 
equations  in  crossflow  plane.  By  using  numerical  schemes  such  as  ADI  {alter- 
nating-direction-implicit) procedures,  they  provide  a means  of  comparing 
solutions  with  those  obtained  from  standard  boundary-layer  procedures  such  as 
the  Box  scheme  used  in  the  present  study  as  well  as  in  our  previous  studies. 

For  three-dimensional  incompressible  flows  in  a curvilinear  orthogonal  coordi- 
nate system,  these  generalized  equations,  which  for  brevity  we  shall  refer  to 
as  PEBL  (parabolic-elliptic  boundary-layer)  equations,  are  given  by 
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Note  that  if  the  underlined  terms  in  (2.1.2)  and  (2.1.3)  are  set  equal  to  zero, 
we  recover  the  usual  three-dimensional  boundary-layer  equations  which  we  shall 
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refer  to  as  the  standard  three-dimensional  boundary- layer  equations.  In  the 
above  equations  h.|  and  h^  are  metric  coefficients  and  are  functions  of  x 
and  0,  that  is, 

h-j  = h^(x,0)  , h^  = h2(x,o)  (2.1.4a) 


The  parameters  and  are  known  as  the  qeodesic  curvatures  of  the 

curves  0 = const,  and  x = const.,  respectively.  They  are  given  by 


h^2  30 


hjfi^  3x 


(2.1.4b) 


At  the  edqe  of  the  boundary  layer,  (2.1.2)  and  (2.1.3)  reduce  to 
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(2.1.5b) 

The  solution  of  the  system  (2.1.1)  - (2.1.4)  requires  both  closure  assump- 
tions for  the  Reynolds  stresses,  and  appropriate  boundary  conditions.  The 
closure  model  will  be  discussed  in  Section  2.4.  Some  of  the  boundary  conditions 
which  need  to  be  applied  depend  upon  whether  the  standard  or  PEBL  equations  are 
to  be  solyed.  In  either  case,  however,  edge  and  wall  boundary  conditions  for 
zero  mass  transfer  are: 


y = 0 u,w  =0  V = 0 (2.1.6a) 

y -►  “>  u -►  Ug(x,0)  w ->  Wg(x,0)  (2.1.6b) 

The  standard  boundary-layer  equations  also  require  starting  solutions  on 
two  intersecting  planes.  One  of  these  is  the  windward  plane  of  synmetry  of  the 
body.  The  other  is  chosen  as  the  line  x - constant  around  the  nose  of  the 
body  where  the  value  of  x is  equal  to  the  value  of  x at  the  stagnation 
point.  On  this  line  the  assumption  is  made  that  the  flow  behaves  as  if  it  were 
along  an  azimuthal  attachment  line.  This  eliminates  the  need  to  trace  stream- 
lines from  the  stagnation  point  around  the  nose  of  the  body.  Since  for  all  the 
bodies  considered  in  this  report  the  stagnation  line  is  extremely  close  to  the 
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nose,  the  effect  of  this  assumption  is  not  felt  after  a few  x steps  downstream 
of  the  starting  solution  due  to  the  parabolic  nature  of  the  boundary- layer 
equations . 

On  the  windward  plane  of  symmetry  w = 0,  but  w.  f 0,  and  all  other 
o-deri vati ves  are  zero.  The  e-momentum  equation,  (2,1.3)  is  singular  on  this 
line,  but  by  differentiating  it  with  respect  to  e , a complete  set  of  equa- 
tions is  obtained.  Thus,  the  windward  symmetry  plane  equations  are; 

Continuity 

+ |y(hih2v)  = 0 (2.1.7) 

^ ^ ^ RJ  '-I  - Vee'^ 
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Here  w^  = 3w/3e,  W0g  = aw^/ae.  These  equations  are  subject  to  the  following 
boundary  conditions 

y * 0 u.Wg  = 0 V = 0 (2.1 .10a) 

y ^ u Ug(x)  Wg  W9g(x)  (2.1.10b) 

The  equations  on  the  azimuthal  attachment  line  are  very  similar  to  the 
symnetry  plane  equations  except  that  now  the  assumption  is  made  that  u = 0, 

Ujj  f 0,  and  all  other  x derivatives  are  zero.  Note  again  that  this  is  just 
an  approximation.  Now  the  x-momentum  equation  is  singular  but  can  be  utilized 
after  differentiation  with  respect  to  x.  This  yields  the  azimuthal  attachment 
line  equations : 

Continuity 

Vx  * Is- • 0 


57  'V*  * ^"a 


X -moment urn 
u 3U 


_ _ . V M - “a  , J_ 

h.|  3X  ay  ITj"  3X  3y 


9 -momentum 

u 1 2 ^*^0 

hNj'  ^ h^  '^6  3y  '^'^9*^2 


5 


r 


X -moment urn 


TO  3U 

1 u2  ^ w X ^ 

fiy  X hg  38  3y 


u wK,  = ^ ^ - u_w  K, 

X 1 h-j  xe  h2  36  xe  e 1 


3 r 

3y  I 


3U 

—A  _ (iTV) 
3y  ' 'x 


z-momentum 


w 3w  3W 

h,  38  3y 


3w„ 

e e 

38 


(2.1.12) 


(2.1.13) 


Here  = 3u/3x,  = 3Ug/3x  and  the  equations  are  subject  to  the  boundary 

conditions 

y = 0 Uj^.w  = 0 V = 0 (2.1.14a) 

y -»•  oo  -*•  Uj^(6)  w -►  w(e)  (2.1.14b) 

The  flow  at  the  stagnation  point  itself  is  the  first  point  calculated  in  the 
solution  procedure.  Here  a combination  of  both  attachment  line  and  symmetry 
plane  procedures  is  used  to  produce  the  equations  for  a general  three-dimensional 
stagnation  point  boundary  layer. 


If  the  PEBL  equations  are  to  be  solved,  rather  than  the  standard  equations, 
it  is  necessary  to  prescribe  conditions  on  two  boundary  surfaces  of  the  azi- 
muthal direction  since  equations  (2.1.2)  and  (2.1.3)  are  elliptic  in  the  e 
coordinate.  On  a body  of  revolution,  two  symmetry  planes  exist  on  the  wind- 
ward and  leeward  generators.  Thus,  at  both  6*0  and  e = it  we  have 

Ug  = 0 ; w = 0 (2.1.15) 


It  is  not  necessary  to  solve  different  symmetry  plane  equations,  as  is  done 
in  boundary-layer  calculations,  to  start  the  calculation  in  the  e-direction. 
The  equations  (2.1.2)  and  (2.1.3)  are  simply  solved  in  the  entire  e-y-plane 
subject  to  the  appropriate  boundary  conditions  stated  above.  An  initial  data 
plane  is  also  required,  as  in  the  standard  case,  for  all  the  calculations 
to  be  made  here  using  the  PEBL  equations.  The  standard  equations  are  used  to 
generate  this  solution,  i.e.,  (2.1.11)  - (2.1.14). 
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2.2  The  Boundary  Layer  on  a Conical  Surface 


If  we  are  specifically  interested  in  boundary  layers  on  conical  surfaces, 
then  we  have  h^j  = 1 and  h2  = x.  It  follows  from  (2.1.4b)  that  = 0 and 
= -1/x.  For  compressible  flow,  the  standard  three-dimensional  boundary- 
layer  equations  equivalent  to  (2,1.1)  - (2.1.3)  are: 

Continuity: 
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Here  e denotes  the  polar  coordinate  in  the  developed  plane;  x the  coordinate 
along  the  generators;  and  y the  coordinate  normal  to  the  surface;  with  w,  u, 

V the  velocities  in  the  e,  x,  and  y directions.  The  boundary  conditions 
for  these  equations  have  already  been  given  as  (2.1.6a)  and  (2.1.6b). 


In  compressible  flow  it  is  also  necessary  to  have  the  energy  equation  to 
close  the  system  of  equations.  In  the  present  notation  It  becomes: 

Energy: 

3H  . A 1 ^ a /u^  + w^  \ 

up*  "V  " 

(2.2.4) 

where  H is  the  stagnation  enthalpy,  and  Pr  the  Prandtl  number.  The 
boundary  conditions  on  the  enthalpy  are 

y ■ 0 H » given  or  aH/ay  « given  (2.2.5a) 

y^-  H^Hg(x,e)  (2.2.5b) 

At  the  windward  stagnation  line,  e * 0,  the  cross-flow  momentum  equa- 
tion is  singular.  Taking  into  account  the  syimietry  conditions  and  differenti- 
ating (2.2.3)  with  respect  to  6,  we  can  write  It  as 
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Noting  that  w •=  0,  (2.2.1),  (2.2.2)  and  (2.2.4)  can  be  simplified  and  written 

as 
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where  w^  = 3w/3e.  The  boundary  conditions  for  this  set  of  equations  is  given 
by  (2.1.10a)  and  (2.1.10b)  plus  (2.2.5a)  and 


H - Hg(x) 


(2.2.10) 


2.3  Geometry  and  Coordinate  System 

Three-dimensional  boundary  layers  develop,  in  general,  on  a nondevelopable 
surface.  This  means  that  the  computations  cannot  be  made  on  a Cartesian  coord- 
inate system;  Instead  a curvilinear  coordinate  system  must  be  used.  There  are 
several  possible  coordinate  systems.  One  possibility  is  a streamline  system. 
Here  the  external  streamlines  form  one  family  of  the  coordinate  lines,  and  the 
equipotential  lines  (orthogonal  everywhere  to  the  streamlines)  form  the  other 
family  of  the  coordinate  lines.  The  use  of  this  system  enables  the  streamwise 
and  crosswise  velocity  profiles  to  be  computed  easily  and  is  one  of  the  most 
popular  coordinate  systems  used  in  integral  methods.  However,  the  determina- 
tion of  the  streamline  coordinate  system  requires  a knowledge  of  the  external 
flow  in  considerable  detail.  In  addition,  the  numerical  computation  of  the 
associated  geometrical  parameters  is  quite  involved.  As  a result,  when  a dif- 
ferential approach  such  as  ours  is  used,  this  coordinate  system  is  not  very 
attractive. 


Another  coordinate  system,  popular  for  bodies  of  revolution,  is  the  one 
used  by  Blottner  and  Ellis  [20].  Here  the  coordinate  system  is  determined  by 
the  intersection  with  the  body  surface  of  parallel  planes  which  pass  through 
an  axis  containing  the  stagnation  ooint.  The  other  coordinate  lines  are 


8 


r 


obtained  numerically  from  the  orthogonality  condition.  When  this  system  is 
employed,  the  computations  can  start  immediately  from  the  stagnation  point, 
pass  through  the  nose  of  the  body  and  march  downstream  to  the  separation  point. 
Again  this  coordinate  system  requires  a lot  of  numerical  work  to  determine  the 
computational  and  the  geometrical  parameters. 

At  small  angles  of  attack  the  stagnation  point  is  close  to  the  nose  of  the 
body.  A very  attractive  coordinate  system  is  the  body  coordinate  system,  in 
which  one  coordinate,  say  e,  is  formed  by  the  lines  x = constant,  and  the 
other  coordinate,  s,  is  the  line  e = constant  along  the  meridional  direc- 
tion, either  measured  along  the  surface  or  measured  along  the  x-axis.  For 
numerical  work  the  latter  is  preferred.  It  is  this  body-oriented  coordinate 
svstem  which  we  shall  use  in  the  oresent  work,  see  Section  5.1  fnr  a snecific 
example. 


2.4  Closure  Assumptions  for  the  Reynolds  Stresses 

The  solution  «f  any  of  the  sets  of  boundary-layer  equations  derived  in 
Section  2.1  requires  closure  assumptions  for  the  Reynolds  shear  stresses, 
-ou'v' , -pw'v' . As  in  our  previous  studies  we  use  the  eddy-viscosity  con- 
cept and  define 


-u 


V-  = 


'1  3y 


(2.4.1) 


Following  the  formulation  of  reference  19,  we  assume  that  e,  = c.  = t 

I 2 m 

and  define  by  two  separate  formulas.  In  the  so-called  inner  region  of 
the  boundary  layer,  is  defined  by 


where 


c = c.  = L 
m 1 


[(tf  * (§fj" 


L^S(y) 


L » 0.4y[l  -exp(-y/A)],  s(y) 

A ” 26  ^ 


u 

T 


(2.4.2) 

(2.4.3a) 

(2.4.3b) 
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In  the  outer  region  is  defined  by 


where 


'■m 


0.0168 


- u^)dy 


= (u~  + wf) 
te  e e 


1/2 


1/2 


(2.4.4) 


(2.4.5) 


where  the  inner  and  outer  regions  are  established  by  the  continuity  of  the 
eddy-viscosity  formula. 
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III.  TRANSFORMATION  OF  THE  GOVERNING  EQUATIONS 


3.1  Transformation  of  the  Standard  Boundary-Laver  Equations 

The  boundary -layer  equations  can  be  solved  when  they  are  expressed  either 
in  physical  coordinates  or  in  transformed  coordinates.  Each  coordinate  system 
has  Its  own  advantages.  In  problems  where  computer  storage  becomes 
important,  the  choice  of  using  transformed  coordinates  becomes  necessary,  as 
well  as  convenient,  since  the  transformed  coordinates  allow  large  steps  to  be 
taken  in  the  x and  e directions.  The  reason  for  this  is  that  the  profiles 
expressed  in  the  transformed  coordinates  do  not  change  as  rapidly  as  they  do 
when  they  are  expressed  in  physical  coordinates.  The  use  of  transformed  coor- 
dinates also  stretches  the  coordinate  normal  to  the  wall  and  takes  out  much  of 
the  variation  in  boundary- layer  thickness  for  laminar  flows.  Since  the 
standard  and  parabolic-elliptic  equations  will  be  solved  by  different  methods, 
and  because  there  are  additional  terms  in  PEBL,  the  transformation  of  each 
set  will  be  considered  separately. 


For  the  standard  equations,  we  first  define  the  transformed  coordinates  by 

.1/2 


X = X 


z = e 


v"i/ 


dy 


(3.1.1) 


and  introduce  a two-component  vector  potential  such  that 


h,u  = ^ 
"2“  3y 


. h,h,v.-(|i.|l)  (3., .2) 


In  addition,  we  define  dimensionless  f and  g related  to  and  4>  by 

'1'  = (Ugvs^)^^^h2  f(x,z,n) 


<P  = ^ g(x,2,n) 

e 


(3.1.3a) 
(3.1.3b) 

Here  Sp  which  denotes  the  arc  length  along  the  x-coordinate,  is  defined  by 

(3.1.4) 


■-/ 


h.|dx 
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Introducing  the  expressions  (3.1.1)  to  (3.1.4)  into  (2.1.2)  and  (2.1.3)  for 
the  standard  equations,  we  get  after  considerable  algebra. 


x-momentum 


(bf")'  + m.ff  -m2(f')^  " -mg(q')^  + m 


m2 


11 


= m 


10 


z*momentum 


, P 2 

(bg”)  + m^fg”  - m3(g')  - m^f'g’  + mggg"  -nig(f’)  + ni^2 

= '"10  (^’  If  ■ I7)  " ^ (9’  If  - 9"  H) 

Here  primes  denote  differentiation  with  respect  to  n and 


f = — 
"e 


g-  =«- 

“e 


b = 1 + e. 


+ 

^m  ^ V 


The  coefficients  to  m^2 

'"1  = 


(3.1.6) 


(3.1.7) 


Si 


Si 


2 = 1177^ 


Si  3Wg 


_ _ e I 6 _ If  c 


■"6  ' Ufi  ^ 


W^  S,  3w  $1  3U 

3z  ^in 


m - 

'"7  - h2  Ug 


n>9  - ^1$! 


(3.1.8) 


■"10  “ ^ 


*"11  ^ '**2  *"5  ^ *"8 


*"12  “ '"4  ■*■  ™3  ■*■  "’9 


To  transform  the  symmetry  plane  equations  (2.1.7)  - (2.1.9),  we  use 
relations  (3.1.1)  and  (3.1.4)  for  the  coordinates,  with  the  dependent  vari- 
ables defined  as 


12 


h,u  . 

2 3y 


N".  ■ 


''iV' -(If*  ♦) 


The  dimensionless  f and  g are  defined  by 

<l>  = (UgvS^)^''^h2f(x,n)  (3. 1.10a) 

* = g(x,n)  (3.1.10b) 

' e 

Using  these  variables  the  symmetry  plane  equations  become 

(bf")  + m^ff"  - m2(f')^  + mgf'q  + = m^Q  (f  1^  - f"  Uj  (3.1.11) 

(bg")  + m^fg"  -m3(g')2  + mggg"  + m^2  = >"10  (^’  !?“-«"  fj) 

(3.1.1?) 


where 


1 


$1 


"'i  ^ ly 


m - -Ve 

"^4  w,.h. 


’ze"l 


3X 


S1K2 


^1 

* I ' — — 

2 Ugh^  az 


mg  = m3 


m.  = 


1‘1  3Wg 


3 Ugh2  3Z 


"10  ■ h 


"’ll  “ *"2 


mi2  - m3  + m^ 


(3.1.13) 


In  the  same  manner,  the  attachment  line  equations  (2.1.11)  - (2.1.13) 
can  be  transformed  by  using 


and 


X = X 


h,u  = 

2 X 3y 


z = e 


h.w  « ^ 
1 ay 


9(z.n) 


(3.1.14) 

(3.1.15) 

(3.1.16a) 

(3.1.16b) 
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With  these  variables,  the  attachment  line  equations  are 

(bf)'  + ff”  -mgfg'  + mggf  + = m^  (g'  ^ - f"  |f) 

(bg")  + fq"  -m3g'^  + m^gq"  + m^^  = "'7  ('>'  “ 9"  |f) 


(3.1.17) 


where 


’3  - hj  3Z 


w 


3U. 


_ , 1,  _e_  f J “'xe 


m “e 


hj  3Z  '"1  “xe  "e' 


(3.1.18) 


'2  “xe 

mil  = 1 + mg 

mi  2 = m3 

The  general  boundary  conditions  which  apply  to  the  three-dimensional  boundary- 
layer  equations  are 

n-0  ^=^w“^w“5w“°  (3.1.19a) 

n = n f = 1.  g’  = 1 (3.1.19b) 


3.2  Transformation  of  the  Conical  Boundary-Layer  Equations 

It  is  well  known  that  for  laminar  flows  with  nonporous  walls  and  wall 
temperatures  independent  of  x,  or  adiabatic  wall  conditions,  the  equations 
(2.2.1)  to  (2,2.4)  have  similarity  solutions  (i.e.,  solutions  Independent  of 
x)  if  we  define  the  similarity  variable  n by 

1/2 

dn  . 1— ody  (3.2.1) 


V'e",*/ 


In  addition,  if  we  introduce  a twc-component  vector  potential  such  that 
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pux 


pwx 


8^ 


3X  X 30 


(3.2,2) 


3y  ' ’'••••  ay 

and  dimensionless  functions  ii>  and  ^ by 

'I'  “ r"  (3.2.3) 

“e 

then  It  can  be  shown  that  the  two  momentum  equations,  (2.2.2)  and  (2.2.3),  and 
the  energy  equation  (2.2.4)  can  be  written  as 

(bf)  + I ff"  + m.,gf  + m2[(g')^  - f g']  = m^^g'  ll^- f |3.j  (3.2.4) 


(bg")  + |.  fg"  + m^gg”  + m3 


+ — - f g' 
P ^ 


"2(5'  if  - 9"  if)  t3-2-5) 

■ ”2(9'  If- If)  (3-9-«) 


+ f fE*  + m^gE' 


where  primes  denote  differentiation  with  respect  to  n. 
The  definitions  of  the  terms  In  (3.2.4)  to  (3.2.6)  are: 


_ U I _ W C k r/i  . 

^ * JT  • 9 - - . E . fp  , b = C(1  + e^) 

C C c 

2 

C ■ • m-i  = 5-  + ST r (p  p ) , 

PgWg  1322  pgPg  00  'Pe^e'  ’ 

m - 1 

"'3  irar 

e 


4'- 

m V 
w. 


” n: 


(3.2.7) 


Here  we  have  used  the  eddy  viscosity  and  turbulent  Prandtl  number  concepts  and 
have  written  the  turbulence  terms  by 


-pu'v'  » pe_ 

p m ay 


-pw'v'  » pe 


3w 
m ay 


(3.2.8a) 

(3.2.8b) 
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Similarly,  equations  (2.2.6)  to  (2.2.9)  can  be  transformed  and  expressed  in 
a form  sii>’,lar  to  those  given  by  (3.2.4)  to  (3.2.6)  by  using  the  transforma- 
tions given  by  (3.2.1)  to  (3.2.3)  with  slight  modification  to  the  following 
parameters , 


pwe  = ly  = (3.2.9) 

(w  ) 

'll  = (PgUgUg)^^^x^^^f(n)  <p  = (PgUgUgX)^'^^  — g(n)  (3.2.10) 

e 

Equations  (2.2.6)  to  (2.2.9)  then  become 

(bf")'  + I ff"  + m^gf”  = 0 (3.2.11) 

(bg")’  + Ifg"  + m3gg"  + (g')^j  + = 0 (3.2.12) 

2 “1  * 

[pF  ('  * ^ -k)  ""]  * I 


Here 


e'e 


The  boundary  conditions  which  are  valid  for  either  (3.2.4)  - (3.2.6)  or 
(3.2.11)  - (3.2.13)  are 


n » 0 f = g = f = q*  = 0 or  given  (3.2.14a) 

n - n,  f - g'  - E =>  1 (3.2.14b) 


The  usual  procedure  that  has  been  followed  in  predicting  the  boundary- 
layer  development  on  cones  at  incidence  is  to  solve  the  system  of  equations 
given  by  (3.2.11)  to  (3.2.13)  for  e - 0,  and  (3.2.4)  to  (3.2.6)  for  e > 0 
for  both  laminar  and  turbulent  flows  with  a mixing-length  or  eddy-viscosity 
formulation.  In  those  calculations  it  is  assumed  that  the  square-root  vari- 
ation of  boundary- layer  thickness  for  laminar  flows  also  applies  for  turbu- 
lent flows. 
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Now,  however,  we  do  not  follow  this  procedure  for  turbulent  flows; 
instead  we  follow  a recent  suggestion  of  Bradshaw  et  al . [4]  and  assume 
that  the  velocity  profiles  u/u^  and  w/w^  are  similar  if  we  define  Y * y/x, 
rather  than  y//x.  It  follows  then  that  for  any  quantity  q 


(3.2.15a) 


(3.2.15b) 


To  have  similarity  = 0;  then  noting  that  aY/Bx  = -y/x^  and 

3Y/3y  = 1/x,  from  (3.2.15),  we  can  write 


(aq/3x)y  (3q/3Y)^(3Y/3x)  ^ 

(3q/3y)x  ' (3q/3Y)^(3Y/3y)  ^ “x  (3.2.16) 

SO  that 


3x  ' X 3y 

With  this  assumption,  it  can  be  shown  that  equations  (3.2.1)  to  (3.2 
be  written  as 

I?  ly  + 2pu  = 0 
2 

w 3u  . 3U  w 3 / 3u  -rrrrr\ 

7 ^ P''  “ iy  iy  ~ ) 

W3W.  ‘3W.  UW_  1^4.3  / 3w  .-7r:7r\ 

^ X — xd?^37(‘^37-p'"'' ; 

X 36  ay  ay  j^Pr  ay  Pr)  ay  \ 5 ) 


- pvTi 


where 


■1 


= £^-u 


Simi 


larly,  the  windward  staqnation-1 ine  equations  can  be  written  as 
pWg  + ly  xpv  + 2pu  » 0 


(3.2.17) 
.4)  can 

(3.2.18) 

(3.2.19) 

(3.2.20) 

(3.2.21) 

(3  .2.22) 

(3.2.23) 

(3.2.24) 
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. wf  UW„  1 .2  , r 3w 

I7  ' ly  [fc  ly  *''('■■  y I7  (r)  “ 


(3  .2.25) 


(3.2.26) 


If  we  now  apply  the  transformation  given  by  (3.2.1)  to  (3.2.3)  with 
pvx  defined  by 


pCx  = -2  i II- 

X X 36 


(3.2.27) 


to  (3. 2. 18)  to  (3.2.21)  and  the  transformation  given  by  (3.2.9),  (3.2.10)  with 
pvx  defined  by 


pvx  = -2 

XX 


(3.2.28) 


to  (3.2.23)  to  (3.2.26),  it  can  be  shown  that  the  resulting  transformed  equa- 
tions are  almost  identical  to  those  given  by  (3.2.4)  to  (3.2.6)  for  the  gen- 
eral case  and  to  those  given  by  (3.2.11)  to  (3.2.13)  for  the  windward 
stagnation  line  case.  The  only  difference  appears  in  the  coefficients  of 
ff",  fg"  and  fE';  instead  of  3/2,  the  new  equations  have  a coefficient 
of  2.  It  should  be  noted,  however,  that  the  resulting  equations  for  turbu- 
lent flows  are  not  strictly  similar  like  the  laminar  flow  equations  because 
they  were  derived  under  the  assumption  of  constant  x.  Thus,  the  equations 
are  valid  for  a fixed  x and  no  scaling  of  solutions  is  allowed. 

3.3  Transformation  of  the  Parabolic-Elliptic  Boundary-Layer  Equations 


Since  the  PEBL  equations  should  reproduce  the  3DBL  results  in 
regions  where  the  cross-flow  gradients  are  not  large,  the  same  argument  can 
be  made  for  using  a similarity  transformation  in  these  equations  too.  In  the 
nondimensional  PEBL  variables  used  in  Appendix  A,  the  transformation  (3.1.1) 
becomes  1 ,0 


/ u \ 

= (r) 


(3.3.1) 


where  u^  * Ug(x,6).  Placing  this  into  the  PEBL  equations  (A. 11)  - (A. 13) 
and  defining  a new  velocity-like  variable,  V as 


18 


V = c 


(3.3.2) 


”x  h2  % 


yields  the  following  equations 
x-momentum 


1 V 2 1 1 1 

uu.  + 7—  wu,  + -T  u - K,uw  + K,w  = Pv  ■*■  D 

£;  h2  z 5 n 1 2 h-j  "^x  5 nn  R ^2  zz 


(3.3.3) 


-momentum 


1 V 21  “e  1/1 

uWf  ■*■  XT-  ww,  + r w - K,uw  +K,'i  = -r-p„+  — w +d’I“tw, 

C h^  z ? n 2 1 9 c on  R I ^2  zz 


— 2 ^ 

^ 2 z 


(3.3.4) 


continuity 


/u 

(— I — 

u 1 ^ 1 

V / 

- K„u  + uy^  J-  (wu^/^)  - K,w  + ]r  V = 
2 e h2  3Z  \ e / 1 c n 


(3.3.5) 


The  appearance  of  the  u factor  in  the  momentum  equations  (3.3.3)  and 
(3.3.4)  is  somewhat  unusual  and  could  be  eliminated  by  normalizing  the 
velocities  with  their  edge  values.  However,  due  to  the  presence  of  the  second 
derivatives  with  respect  to  z now  present  in  the  equation,  this  would  intro- 
duce more  terms  than  are  normally  encountered  (such  as  a Ue^^  term)  and 
it  was  deemed  advisable  to  leave  the  equations  as  presented  above.  This  led 
to  a computational  problem,  to  be  discussed  in  Section  4.2,  which  was  solved 
by  a new  technique  also  described  in  Section  4,2. 

The  boundary  conditions  to  be  imposed  on  this  transformed  set  of  equa- 


tions  are 

n = 0 

c 

II 

It 

o 

V = 0 

(3.3.6a) 

n » n 

CO 

U - Ug(t,z) 

w ►Wg(c,z) 

(3.3.6b) 

Z « 0,n 

o 

II 

II 

o 

(3.3.6:) 

IV,  NUMERICAL  METHODS 


4.1  Numerical  Formulation  of  the  Standard  Boundary-Layer  Equations 

We  use  the  Box  method  to  solve  the  governing  3DBL  equations.  This 
is  a two-point  finite-difference  method  that  has  successfully  been  applied 
to  two-dimensional  flows  by  Keller  and  Cebeci  and  to  three-dimensional  flows  by 
Cebeci.  A detailed  description  of  the  method  is  presented  in  [1]  and  [21]. 

For  this  reason  only  a brief  description  of  it  will  be  presented  here. 

One  of  the  basic  ideas  of  the  Box  method  is  to  write  the  governing  sys- 
tem of  equations  in  the  form  of  a first-order  system.  Thus,  in  our  case,  the 
first  derivatives  of  f and  q with  respect  to  n are  introduced  as  new 
unknown  functions.  With  the  resulting  first-order  system  and  an  arbitrary 
rectangular  net,  we  use  simple  centered  difference  quotients  and  averages  at 
the  midpoints  of  net  rectangles  and  net  segments  to  get  second-order  accurate 
finite-difference  equations.  Then  nonlinear  difference  equations  are  linear- 
ized by  using  Newton's  method  and  the  resulting  linear  system  is  solved  by 
the  block-elimination  method  discussed  by  Isaacson  and  Keller  [22]. 

In  our  present  method  we  solve  the  two  momentum  equations  simultaneously. 
Essentially  the  stagnation-line  equations  and  the  symmetry  plane  equations  are 
two-dimensional  flows  in  the  sense  that  these  equations  have  only  two  independ- 
ent variables,  (x,a)  or  (z,ri).  On  the  other  hand,  the  two  momentum  equa- 
tions (3.1.5)  and  (3.1.6)  are  three-dimensional  flow  equations  for  obvious 
reasons.  The  solution  of  the  two-dimensional  flow  equations  is  discussed  in 
considerable  length  in  references  [18],  [21  ],  for  this  reason  we  shall 
only  discuss  the  solution  of  three-dimensional -flow  eouations,  namely,  (3.1.5), 
(3.1.6)  and  their  boundary  conditions,  (3.1.19). 

With  the  introduction  of  new  independent  variables  u(x,z,n).  v(x,z,n), 
w(x,z,n)  and  t(x,z,n).  the  equations  given  by  (3,1.5)  and  (3.1.6)  can  be 
written  as 


f 

= u 

(4.1.1a) 

u' 

= V 

(4.1.1b) 

g' 

= w 

(4.1.1c) 

w' 

= t 

(4.1. Id) 
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(bv)'  * m,fv  - - mjUW  ♦ m^vg  - |i  _ , ^ j 

(bt)'  + m^ft  -m^uw  - m3W^  + m^qt  - tiigU^  + = "'lO  1?"^  fx') 

We  next  consider  the  net  cube  shown  in  Figure  1 and  introduce  the  net 
points  by 


= 0 

"n  = "n-l  * 

n = 1, 

2. 

....  N 

= 0 

■ 'i-l  ♦ '■i 

i = 1. 

2, 

....  I 

(4.1.2) 

= 0 

"j  ■ "j-1  * "j 

j = 1. 

2, 

....  J 

The  difference  equations  which  are  to  approximate  (4.1.1a)  to  (4.1. Id) 

are  obtained  by  averaging  about  the  midpoint  (x  , z.,  n.  wo) 

n’  1 j-1/2' 


Figure  1.  Net  cube  for  the  difference  equations  for  three-dimensional  flows. 
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(4.1.3a) 


(4.1.3b) 

(4.1.3c) 

(4.1.3d) 


The  difference  equations  which  are  to  approximate  (4.1.1e,f)  are  rather 
lengthy.  To  illustrate  the  differencing  of  equations  similar  to  (4.1.1e,f), 
we  consider  the  following  model  equation: 

V = "’lo"  It  ^ V If 

The  difference  equations  for  this  equation  are 


* ("7)m/2-j.1/2(^^-)(«''-5) 


where,  for  example. 


4- ('S'' 

-t'-S''  4“3‘’’’  *“3:!  4“S-1’’> 
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The  boundary  conditions  for  the  system  of  equations  given  by  (4.1.1)  at 

X = x_  and  at  z = Zj  are: 
n 1 


n,i 

-n,i  n n,i  . n-i  - n,i  » n,i  , n,i  I \ 

Fq-  =0,  Pq*  =0.  Uq  =0.  Wq*  = 0.  Uj  =1.  Wj 


If  we  assume  u? 

J 


pn,i-l  ,n,i-l 


f Uj  * J * * 3nd  * 

to  be  known  for  O^j  ^J,  then  the  differ- 
ence equations  (4.1.3)  and  the  difference  equations  for  (4.1.1e,f)  along  the 
(4.1.6)  yield  an  implicit  nonlinear  algebraic  system  of  6J  + 6 equations  in 
as  many  unknowns  (f^,  Uj,  Vj,  g^,  w^,  tj).  We  solve  this  nonlinear  system 
by  means  of  Newton's  method.  The  resulting  linearized  system  is  then  solved 
by  using  the  block  elimination  method  discussed  by  Isaacson  and  Keller  [22]. 


n-l,i-l  n-l,i-l  n-l,i-l 


(«.1.6) 

n-l,i-l  ,.n-l,i-l 


n,i-l  *n,i-i 


'J 


n-l,i 


t'J 

n-1  ,i 


4.2  Numerical  Formulation  of  the  Parabolic-Elliptic  Boundary-Layer  Equations 

The  PEBL  equation  set  (3.3.3)  - (3.3.5)  is  elliptic  in  the  cross-flow 
Z-T1  plane,  but  still  allows  a parabolic  march  in  the  main  flow,  direction. 
An  implicit  numerical  procedure  was  chosen  to  integrate  these  governing  equa- 
tions for  a number  of  reasons.  The  success  of  implicit  methods  on  the  boundary- 
layer  equations  and  parabolic-elliptic  Navier-Stokes  equations  [23]  implies 
that  they  should  be  efficient  for  the  very  similar  PEBL  equations.  It  is 
expected  that  solutions  will  be  required  at  arbitrary  points  along  the  body; 
and  so  for  convenience  in  the  computation,  to  eliminate  the  need  to  check  the 
step  size  restrictions  of  explicit  methods,  unconditionally  stable  methods 
(which  are  consistent  with  the  original  partial -differential  equations  of  the 
problem)  are  necessary.  This  leads  to  the  consideration  of  implicit  methods. 

The  particular  method  used  in  this  study  is  the  alternating-direction 
implicit  (ADI)  method,  used  in  the  form  originally  proposed  by  Peaceman  and 
Rachford  [24].  The  ADI  method  is  ideally  suited  for  the  solution  of  (3.3.3)- 
(3.3.4).  There  is  no  stability  restriction  on  the  step  size,  and  hence 
arbitrary  t-steps  are  permitted.  The  method  has  second-order  truncation 
error  in  its  marching  variation,  which  is  also  a requirement  for  the  type  of 
flow  envisioned  here,  since  the  L-history  of  the  flow  must  be  traced  accurately 
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at  each  step.  Finally,  the  method  does  not  require  the  inversion  of  a sparse 
banded  matrix,  as  a fully  implicit  or  Crank-Nicolson  scheme  would.  Simple 
tridiagonal  coefficient  matrices  are  generated  at  each  step  which  require  much 
less  storage  and  time  for  their  inversion  when  compared  with  sparse  matrices. 


The  ADI  method  is  a two-step,  second-order  accurate  procedure  centered 
about  the  midpoint  of  the  c step  from  iAC  to  (i  + 2)ac,  i.e.,  at 
(i+l)A5,  for  uniform  az  and  Ari  grids.  An  example  of  its  use  is  shown 
on  the  following  equation  modelling  either  (3.3.3)  or  (3,3.4). 


"l“x 


"2“y 


"4  * “zz  * “ 


yy 


(4.2.1) 


The  two  ADI  steps  are  written 


AX  ^^l^j,k  ^“j,k  “j.k^  ^ ^^2^j,k  Vj.k  ^ ^'^3^j!k  ^z'^j.k 


i+1  , ..i 


J+l 


(4.2.2a) 


" ^ Vi.k 


(4.2.2b) 


= <"4>j!l  * 


+2 

k 


where  the  m ^3  are  nonlinear  coefficients,  is  a source  term, 

and  the  standard  central  differencing,  including  the  option  of  nonuniform 
spacing  has  been  used,  i.e., 

, (44.)S.k*l  - 

■ ' ■'  rAZ^)(AZ_)l(AZ^)  + (AZJJ 

(4.2.3a) 


Vj.k 


2 i tl44.)l.{.|,.i  -[(4Z.)  X (4Z.)]u'  * (4Z.)uj_;,.,  ) 

'^Z  j,k  I Caz  )(az  )L{az  ) + (az  )J  I 

The  pressure  gradient  portion  of  N^,  since  it  is  a known,  prescribed 
function,  given  by  the  external  flow  at  all  the  marching  steps,  is  computed 
from 


Px)^*’«1CPx)^^Px)’*'] 


(4.2.4) 
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and  not  simply  evaluated  at  the  (*)  station,  which  is  only  a temporary 
integration  location,  and  does  not  correspond  to  the  physical  midpoint  (i  + 1). 
The  remaining  are  given  by  a second-order  extrapolation  (of  consistent 
accuracy  with  the  integration  procedure)  from  previously  computed  x-stations. 
Originally,  the  standard  extrapolations  were  used.  These  are  derivable  from 
Taylor’s  series,  i.e.. 


{4.2.5a) 


(4.2.5b) 


depending  upon  where  the  variables  are  stored.  However,  since  the  values  of 
the  dependent  variables  used  in  the  PEBL  equations,  (3.3.3)-(3.3.5)  were  not 
normalized  by  their  edge  values,  this  procedure  yields  spurious  results  near 
the  edge  of  the  boundary  layer  and  in  the  inviscid  portion  of  the  flow.  A 
simple  correction  for  this  difficulty  has  been  found  [25],  and  incorporated 
into  the  PEBL -ADI  integration  procedure.  It  consists  of  incorporating  the 
known  information  about  the  edge  velocities  into  the  linearization  so  that  it 
correctly  predicts  the  midpoint  values  at  (i  + 1).  This  is  accomplished  by 


C = 


[(Qe  Qe^)/2] 

[(1  + AX^/2AX_)Qg  - (Ax^/2iix_)Qg’^ 


(4.2.6) 


where  the  Q in  (4.2.6)  is  the  usual  value  given  in  (4.2.5a). 


In  this  way,  the  two  momentum  equations  are  decoupled  from  each  other 
and  the  continuity  equation  (3.3.5),  and  the  integration  can  proceed  without 
iteration.  After  the  complete  two-step  ADI  procedure  is  completed,  the  solu- 
tion obtained  is  second-order  accurate,  centered  at  (i  + 1)ac.  It  is  at  this 
^•station  that  the  continuity  equation  is  evaluated.  The  unknown  velocity 
V(5,2,n)  can  be  computed  at  the  given  ^-station  successively  out  from  the  wall, 
where  the  boundary  condition  V(l,z,0)  = 0 holds,  to  the  edge  of  the  computa- 
tional domain,  and  around  the  cone  from  the  windward  to  leeward  meridian. 

This  is  accomplished  by  differencing  equation  (3.3.5)  as  follows 
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where 


* hj«z'^j-i/2.k  -I'l  " h7;;TTT/'^j-i/2.k 


i+l 


.iH 


Vj-l/2,k  -\h  * i:T^fj-l/2.k 


.1+1’ 


- K,  + ^ 


.1+1 


(4.2.7) 


6 

%''j-l/2,k 


„i+l  _ .1+1 

J.k  ^j-l>k 


An 


j 


T .1+1 

^C“j-l/2.k 


1+2  _ 1 
^.1-1/2.k  “j-1/2.k 

2 At 


6^  is  defined  by  (4.2.3a)  and  all  the  Indicated  differentiations  In  (3.3.5) 
have  been  performed. 


The  differencing  (4.2.2)  generates  tridlagonal  matrices  for  each  unknown 
at  both  values  of  the  ADI  step* which  are  easily  solved  by  the  familiar 
algorithm.  The  only  alteration  from  the  straightforward  application  of  the 
integration  takes  place  on  the  z-implicit  step  for  the  u-velocity.  Here,  the 
boundary  conditions  to  be  applied  are  zero  gradient  at  both  boundaries,  see 
(3.3.6c).  Rather  than  differencing  the  equation  on  the  boundaries  and  using 
fictitious  nodes  to  fulfill  boundary  conditions,  the  integration  domain  was 
limited  to  include  all  but  the  boundary  points.  The  boundary  condition  was 
then  incorporated  by  using  a second-order  accurate  one-sided  difference  to 
define  the  values  of  the  first  and  last  points  in  terms  of  its  nearest 
neighbors  as  follows: 


^^2^3, k 


u = ^^"2  * ""3^  “2.k 

[(aZ2  + AZj)  -AZ2] 

2 2 

..  ^ ^^JMAX-1^  “jMAX-l,k  “ ^^JMAx“jMAX-2,k 

'jMAX,k  r..  ’2  ^ 

JMAX  ^^JMAX-r  ^^JMAxJ 


(4.2.8) 


The  procedure  does  not  destroy  the  tridiagonal  nature  of  the  coefficient 
matrix,  and  has  been  used  with  success  previously  [23]. 
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V.  RESULTS 


5.1  Body  of  Revolution 

The  test  case  chosen  as  a standard  for  all  the  body  of  revolution  compar- 
isions  to  be  given  below  is  the  flow  past  a prolate  spheroid,  as  shown  in 
figure  2,  whose  axis  ratio,  t,  equals  1/4. 

The  solution  to  the  three-dimensional  boundary -layer  equations,  written 
in  a general  orthogonal  body-oriented  coordinate  system  requires  the  prescrip- 
tion of  the  body  geometry,  and  derivable  from  that,  the  external  potential  flow. 

In  order  to  assess  the  accuracy  of  the  boundary-layer  solution,  an  analytically- 
prescribed  body  with  reasonably  simple  closed  form  solutions  for  the  potential 
flow  was  used.  This  avoids  the  problems  caused  by  trying  to  spline  fit  and 
smooth  numerically  computed  surfaces  and  velocity  components.  It  must  be 
emphasized  that  this  simple  body  is  not  necessary  for  the  successful  use  of 
the  computer  program,  and  general  geometry  packages  are  available  [1  ].  Thus, 
the  calculations  were  performed  over  a simple  body  of  revolution,  namely,  the 
prolate  spheroid.  This  configuration  was  previously  considered  by  Wang  [12]  - [14] 
who  has  documented  the  resulting  flow  field  at  various  angles  of  attack.  This 
body  has  extensive  regions  of  cross-flow  reversal  within  a still  unseparated 
region  and,  in  addition,  has  an  unusual  separation  pattern  at  moderate  angles 
of  attack  which  makes  it  an  excellent  test  case  on  which  to  perform  three- 
dimensional  boundary-layer  studies.  It  is  felt  that  the  extension  to  compres- 
sible flow  will  not  be  more  than  a minor  alteration  to  the  procedures  developed 
here.  Similarly,  the  inclusion  of  turbulence  can  be  easily  accomplisned  for 
three-dimensional  boundary  layers  by  the  means  of  the  eddy-viscosity  model 
described  in  Section  2.4,  and  used  in  the  following  section  for  conical  flows. 

What  is  essential  is  the  verification  of  the  methods  to  be  used  for  such 
computations.  Hence  this  study  dealt  with  incompressible  flow  in  order  to 
determine  the  features  of  the  two  numerical  techniques  proposed. 

For  the  prolate  spheroid  depicted  in  figure  2 we  can  write  the  equation 
of  the  surface 


Fiqure  2,  Geometry  of  the  prolate  soheroin 


and  the  incremental  distance  along  the  surface 

dsj  = dr^  + dx^ 


and  using  eq.  (5.1,1)  gives 


C^)  - H 


x^{t^  - 1 


(1  - x") 


But,  from  fiqure  2,  we  have 


Hence 


(5.1.2) 


COSB  = 


Q-- x-I  I 

1 + x^{t^  - 1)J 


sins  = 


The  metric  coefficients  are  determined  from 
ds^  = iV|dx 

thus,  from  fiqure  2 and  eq.  (5.1.1), 


[l+x'^'t  -1)]  (5.1.3) 


ds^  = l^de 


ds^  = rde  = t(l  - x^)^^^do 


and  from  eq.  (5.1.2) 
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Consequently, 


1 x^(t^  - 1) 

(1  - x^)  . 


and  the  curvatures  are 


h2=  t(1  (5.1.4) 


[1  + x^(t^  -1)]  (1  -x^) 


= 0 


(5.1.5) 


These  geometrical  relationships  are  used  in  the  governing  equations  of 
boundary-layer  flow  and  also  in  the  determination  of  the  potential.  From 
Faulkner  et  al . [26],the  velocity  components  on  the  surface  of  a prolate 
spheroid  are  given  by 

y — = VQ(t)  cosa  C0S6  - VgQ(t)  sina  sine  cose  (5.1.6) 


5^'^®  (5.1.7) 


where  the  various  angles  a,  6 and  e are  depicted  in  figure  2,  and  the 
functions  Vq  and  '^90  are  only  functions  of  t,  the  thickness  of  the  body, 
given  by  [20] 


Vo(t) 


Vgo(t) 


^ - 1/2  In  [(1  + -^1^  t^)] 

2Vo(t) 

“ 2vTtT-l 

0 


(5.1.3) 

(5.1.9) 


Equations  (5.1.4)  - (5.1.9)  are  all  the  input  necessary  to  the  boundary-1  aver  equa 
tions  for  the  flow  past  a prolate  spheroid. 


The  initial  calculation  made  was  simply  the  prolate  spheroid  at  zero 
degrees  incidence,  i.e.,  axisymmetric  flow.  While  this  is  a trivial  case 
for  the  standard  3DBL  procedure  (it  calculates  one  attachment-line  flow,  and 
the  result  is  valid  at  every  angular  position  around  the  body),  it  has  to  be 
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dencnstrated  that  the  PEBL  procedure,  which  by  its  very  nature  must  calculate 
all  the  points  around  the  body,  generates  axi symmetric  flow.  Consequently, 
the  PEBL  procedure  was  used  for  the  axisymmetric  case  using  only  five  points 
around  the  body,  i.e.,  AO  “ 0.78539  (45°),  a very  crude  grid.  The  grid 

spacings  across  the  boundary  layer  and  down  the  body  were  identical  to  those 
used  previously  for  3DBL,  i.e..  An  = 0.3  and  ax  varying  up  to  0.1  for 
X > -0.5.  As  in  all  cases  to  be  discussed,  the  initial  data  planes  necessary 
for  the  extrapolation  in  (4.2.5)  were  obtained  from  the  standard  boundary-layer 
calculation. 

The  results  of  this  calculation  showed  that  the  initially  axisymmetric 

boundary  layer  remained  that  way  until  the  separation  point  was  reached.  The 

computed  values  of  w,  which  are  exactly  zero  for  axisymmetric  flow,  were 

-8 

all  of  the  order  of  10  (below  the  round-off  error  of  the  single  precision 
arithmetic  of  the  IBM  machine  used).  The  differences  in  the  u velocity 
around  the  body  occurred  in  the  fifth  decimal  place,  again  close  to  round-off 
error. 

The  plot  shown  in  figure  3,  gives  a comparison  of  f",  the  wall 
shear  parameter,  obtained  by  both  methods.  The  curves  are  effectively  indis- 
tinguishable until  the  effects  of  the  adverse  pressure  gradient  at  the  rear  of 
the  body  are  felt.  Here,  one  would  expect  the  lack  of  iteration  in  the  PEBL 
procedure  to  generate  less  accurate  results  than  the  3DBL  procedure  which  uses 
a Newton  iteration.  Nevertheless,  the  results  are  fairly  good,  and  the  axi- 
symmetry  is  indeed  maintained  by  the  PEBL  procedure  only  by  the  imposition  of 
symmetry  conditions,  c.f.  (3.2.6c),  at  the  windward  and  leeward  meridians. 

The  3DBL  results  were  next  compared  with  the  published  results  of  Wang 
for  the  axisymmetric  case  [12].  The  skin  friction  (in  Wang's  notation)  down 
the  body  computed  by  both  methods  is  shown  in  figure  4.  There  is  a notice- 
able difference  between  the  two  calculations  with  separation  predicted  at 
X 2t  0.675  and  x 0.8  by  the  3DBL  or  PEBL  calculation  and  Wang,  respectively. 
An  independent  calculation  of  this  same  problem  made  by  Chang  [27],  substanti- 
ates our  results. 

The  next  test  case  considered  was  chosen  to  generate  just  a small  amount 
of  crossflow.  The  incidence  of  the  flow  was  taken  to  be  2°.  The  general 
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structure  of  the  flow  near  the  body  is  shown  in  fig.  5.  The  ordinate  repre- 
sents the  azimuthal  position  around  the  body,  and  the  abscissa  is  the  x 
coordinate  measured  along  the  major  axis.  The  body  surface  is  thus  the  plane 
of  the  figure.  Results  from  two  separate  calculations,  in  addition  to  some 
points  extrapolated  from  Wang's  study  [12]  are  shown  on  this  figure.  The 

triangles  represent  calculations  made  using  the  3L'BL  method.  To  the  left  of 

the  triangles,  both  components  of  the  shear  stress,  t -v  &u/3y  and 

"v  9w/3y,  are  positive.  The  open  triangles  represent  the  grid  points  at 
which  the  crossflow  shear,  reverses  sign,  as  calculated  by  3DBL.  Hence 

the  region  covered  by  the  triangles  represents  a region  of  reverse  crossflow 
in  the  boundary- layer  flow  past  the  prolate  spheroid.  Computed  separation 
takes  place  at  the  points  given  by  the  solid  triangles.  Beyond  130°  the  grid 
was  too  coarse  in  x to  detect  any  separation  for  x > 0.4.  Note  that  the 
locus  of  separation  points  in  some  sense  parallels  the  line  of  reverse  cross- 
flow 


This  configuration  was  next  calculated  using  the  PEBL  procedure.  Again 
on  Fig.  5,  the  open  circles  represent  the  first  locations  at  each  x-station 
where  < 0.  The  locus  of  these  points  has  been  approximated  by  the  curved 
line.  It  is  obvious  that  it  coincides  with  the  prediction  of  the  3DBL  calcu- 
lation for  the  onset  of  reverse  crossflow.  There  is,  however,  a large 
discrepancy  in  the  predicted  separation  point.  The  PEBL  procedure  first 
senses  separation  along  the  x = 0.65  line  almost  on  the  entire  lee  surface, 
130°  16;^  180°,  see  straight  line  on  Fig.  5.  This  includes  a region  wiiich 
was  not  accessible  to  the  3DBL  calculation  due  to  its  predicted  separation. 

It  should  be  noted  that  the  region  x > 0.65  computed  by  the  3DBL  method  for 
0 ^40°  was  not  accessible  to  PEBL.  The  discussion  of  the  inability  of  the 
3DBL  to  compute  the  entire  lee-side  region  will  be  postponed  until  the  higher 
angle  of  attack  case  is  considered.  The  region  computed  by  3DBL  which  was  not 
computed  by  PEBL  is  not  calculable  by  the  standard  PEBL  procedure  set  forth 
here.  Since  all  the  e stations  are  coupled  together  during  the  calculation, 
when  one  station  encounters  separation  and  cannot  advance  into  the  reverse 
flow,  it  automatically  precludes  the  further  advance  of  all  the  other  stations, 
i.e.,  only  the  most  upstream  separation  point  will  be  predicted.  Thus,  as 
will  be  seen  again  later,  one  point  of  < 0 stops  the  calculation.  This 
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Fiqure  5.  Calculated  surface  shear  distributions  for  prolate  spheroid  at  2° 
incidence;  see  text  for  symbol  key. 
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is  not  the  case  for  the  3DBL  where,  as  long  as  the  separation  line  sweeps 
upstream  in  x in  the  direction  of  the  e integration,  the  calculation  can 
continue. 

To  test  the  PEBL  calculations  on  their  behavior  approaching  separation, 

2 

Fig.  6 shows  a plot  of  against  x.  This  should  give  a linear  curve 
approaching  separation.  Stations  at  160°  and  170°  are  shown  as  the  most 
relevant  since  there  the  u velocity  flow  will  most  nearly  be  perpendicular 
to  the  predicted  separation  line.  The  figure  shows  that,  at  least  at  these 
two  stations,  the  calculation  reproduces  a square-root  singularity  at  sepa- 
ration. A further  check  on  the  calculations  can  be  made  by  using  values  of 
the  shear  on  the  two  symmetry  planes  for  a = 2°  interpolated  from  Wang's 
calculations.  The  open  squares  in  Fig.  5 show  the  location  of  crossflow 
reversal  and  the  solid  squares  indicate  separation.  At  6 » 0°,  the  3DBL 
agrees  with  Wang's  results,  and  at  e = 180°,  the  PEBL  agreement  is  also 
good  despite  the  lack  of  iteration. 

Comparisons  of  the  streamwise  variation  of  the  two  components  of  shear 
for  both  3DBL  and  PEBL  calculations  are  given  in  Figs.  7 and  8.  It  is  diffi- 
cult to  make  any  clear  judgements  between  the  methods.  When  the  crossflow 
is  unidirectional,  there  is  little  variation  in  f"  until  crossflow  reversal 
occurs.  There  the  variation,  and  even  the  trend  in  some  cases,  increases. 

For  the  more  sensitive  crossflow  g"  as  the  adverse  pressure  gradient 
increases  for  x > 0,  the  disagreement  increases. 

One  further  point  can  be  made.  The  separation  and  crossflow  reversal 
patterns  computed  for  this  low  angle  of  attack  case,  a = 2° , fit  within  the 
originally  conceived  structure  set  forth  by  Wang  in  his  early  paper  [l2]. 

When  the  angle  of  attack  is  increased  from  2°  to  6°,  new  features  of  the 
flow  emerge,  and  a clearer  understanding  of  the  reasons  for  some  of  the 
behavior  computed  for  the  a = 2°  case  can  be  obtained.  The  first  general 
result  to  consider  is  again  the  distribution  of  the  surface  shear  over  the 
entire  body;  see  Fig.  9.  As  for  the  2°  case,  this  figure  is  a composite  of 
two  separate  calculations,  PEBL  and  30BL,  and  also  includes  data  from  Wang 
who  gave  a detailed  map  of  the  shear  for  this  6°  incidence  case  [14].  The 
line  of  crossflow  reversal,  as  computed  by  Wang,  connects  the  open  squares; 
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Figure  7.  Streatrwise  surface  shear  variation  on  prolate  spheroid  at  2' 
incidence  for  three  azimuthal  locations. 


Figure  8.  Cross-stream  surface  shear  variation  on  prolate  spftieroid  at  2° 
incidence  for  three  azimuthal  locations 
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surface  shear  distributions  for  prolate  spheroid  at  6 
see  text  for  symbol  key. 
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the  oddly  shaped  line  of  separation  connects  the  solid  squares.  The  cross- 
flow  reversal  (open  triangles)  computed  by  30BL  approximates  the  solid  line 
fairly  well;  much  as  was  the  case  for  a = 2°.  The  line  of  separation  as 
computed  by  3DBL,  the  closed  triangles,  now  unquestionably  follows  the  cross- 
flow  reversal  line,  as  was  suspected  in  the  2°  case  also,  compare  Fig.  5. 

To  determine  whether  this  pattern  is  truly  the  case,  we  have  the  PEBL  results 
and  the  results  of  Wang  as  a comparison.  Referring  to  Fig.  9,  both  other 
methods  agree  on  the  general  shape  and  location  of  the  crossflow  reversal 
line,  but  the  discrepancy  in  the  true  separation  line  is  considerable  once  a 
significant  region  of  reverse  crossflow  has  been  traversed.  Thus,  it  appears 
that  the  flow  computed  using  the  3DBL  procedures  described  in  Section  4.1  gives 
results  which  are  highly  dependent  on  the  position  of  crossflow  reversal. 

This  is  not  unexpected  as  it  has  been  shown  by  many  people,  e.g.,  Krause  [28], 
that  the  z-differencing  shown  in  Fig.  1 and  used  in  (3.1.5),  (3.1.6)  is  unstable 
for  any  negative  crossflow.  This  did  not  obviously  manifest  itself  in  the 
lower  angle  of  attack  case,  but  as  will  be  seen  in  the  following  results  to 
be  presented,  the  instability  is  clear.  The  PEBL  calculation  which  proceeds 
without  difficulty  through  the  reverse  crossflow  finally  computed  separated 
flow  at  one  point  x = 0.375,  e = 115°,  and  stops.  As  was  the  case  for  the 
a = 2°  calculation,  this  point  is  upstream  of  the  computed  separation  line 
of  Wang,  as  has  been  the  case  in  all  the  previous  calculations.  Although 
separation  did  not  occur  until  further  downstream,  the  PEBL  procedure  did 
calculate  unrealistic  behavior  in  the  region  120°  i e ^ 140°  for  x > 3.0 
which  tends  to  confirm  the  upstream  intrusion  of  the  line  of  separation  com- 
puted by  Wang  and  indicated  on  Fig.  9.  Much  of  this  can  be  made  more  clear 
from  the  following  series  of  figures. 

The  variation  of  the  two  components  of  wall  shear  are  shown  in  figures 
10  and  11.  In  all  cases  the  agreement  between  3DBL  and  PEBL  is  good  along 
the  entire  body  until  the  onset  of  crossflow  reversal.  Beyond  this  point, 
the  shear  values  deviate  considerably  with  the  3DBL  computed  values  showing 
an  abrupt  break  in  the  previous  smooth  trend.  The  PEBL  calculations  continue 
with  no  difficulty  until  the  adverse  pressure  gradient  makes  the  extrapola- 
tion process  somewhat  questionable.  Close  to  separation,  the  computed  values 
of  f"  between  120°  i 9 i 140°  show  a minimum  and  then  begin  to  rise.  This 
can  be  seen  in  the  9 - 135°  plot  in  Fig.  10  and  in  Fig.  12.  This  behavior 
is  certainly  unrealistic  and  indicates  a breakdown  of  the  PEBL  computation 
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Figure  11.  Cross-stream  surface  shear  variation  on  prolate  spheroid  at  6° 
incidence  for  four  azimuthal  locations. 
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in  the  odd  interior  region  of  the  flow  where  Wang  predicts  separation.  A 
negative  u velocity  is  sensed  by  the  program  at  e = 115°,  x = 0.375,  and 
the  calculation  stops.  However,  prior  to  this  unrealistic  behavior,  a plot 

p 

of  vs  X shows  that  the  calculated  values  are  approaching  separation  in 
the  correct  way  (with  a square-root  singularity)  along  a broad  front  of  e 
just  on  the  lee  side  of  the  body,  see  Fig.  12. 

A different  manner  of  looking  at  the  results  which  brings  out  other 
characteristics  of  the  computed  flow  is  seen  in  Figs.  13a  - 13c.  Here  the 
azimuthal  variation  of  the  shear  is  given  at  various  stations  down  the  body. 

The  first  figure,  13a,  shows  the  shear  at  a point  on  the  body  just  prior  to 
any  crossflow  separation,  x = -0.3.  The  agreement  between  both  calculations 
is  excellent  as  it  should  be  since  there  is  no  reason  for  either  to  be  suspect. 
Further  downstream,  at  x = 0,  the  calculations  indicate  the  final  third  of 
the  body  on  the  lee  surface  to  contain  reversed  crossflow  as  seen  in  Fig.  9. 
Figure  13b  indicates  good  agreement  between  both  3DBL  and  PEBL  calculations 
until  this  region  but  shows  that  although  the  PEBL  calculation  is  unaffected, 
the  30BL  results  begin  to  oscillate  until  the  computer  program  senses  either 
separation  or  nonconverqence  of  the  iterated  solution.  The  final  x location 
displayed  is  at  x = 0.3,  see  Fig.  13c.  Here  the  crossflow  reversal  exists 
on  a greater  portion  of  the  body  surface,  and  the  3DBL  calculation  proceeds 
less  far  around  the  body  while  the  PEBL  calculation  continued  to  the  lee  surface 
without  difficulty  (not  shown  on  the  scale  given).  No  oscillation  occurs  as  the 
program  stops  computing  almost  immediately  after  the  crossflow  changes  sign. 

The  abrupt  drop  at  e ■ 90  is  attributed  to  the  fact  that  the  Ae  grid  changed 
at  the  same  point  where  the  velocity  also  changed  sign. 

The  behavior  of  the  flow  at  these  three  x-stations  can  also  be  seen  on 
the  velocity  profiles  just  off  the  windward  and  leeward  meridians.  These  are 
shown  in  Figs.  14a  - 14c.  At  x = -0.3,  a basically  standard  set  of  velocity 
profiles  for  u/u^  and  w/w^  is  generated.  The  boundary  layer  is  thicker 
near  the  lee  side,  and  no  unusual  behavior  is  noticed  except  in  the  w pro- 
file at  175°,  which  has  a slight  inflection  due  to  the  adverse  azimuthal 
pressure  gradient.  By  x » 0 the  crossflow  on  the  lee  side  c./ntains  a small 
layer  of  reversed  flow.  The  streamwise  flow  on  the  lee  side  has  become  less 
full,  and  the  thickness  of  the  boundary  layer  here  has  increased.  On  the 
windward  side,  the  shape  and  thickness  of  both  u and  w profiles  have 
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Figure  13.  Continued  (b)  x « 0, 
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changed  hardly  at  all.  By  the  time  x = 0.3  has  been  reached,  the  reverse 
flow  layer  on  the  lee  side  occupies  a considerable  portion  of  the  boundary 
layer  which  has  again  increased  significantly  in  thickness.  There  is  some 
slight  difference  in  the  u profiles  between  x = 0 and  x » 0.3,  but  the 
major  alteration  has  been  in  w.  Again,  on  the  windward  plane,  except  for  a 
slight  increase  in  thickness,  both  profiles  are  relatively  unchanged. 

One  of  the  conclusions  reached  by  Wang  from  his  calculations,  was  that 
there  was  no  evidence  of  vortical  flow  on  the  lee  side  of  the  body  as  the 
flow  approached  separation.  It  was  thought  that  this  might  have  been  due  to 
the  use  of  the  boundary-layer  equations  themselves  and  not  truly  representa- 
tive of  the  flow  field.  The  PEBL  procedure  has  already  been  demonstrated  as 
being  capable  of  capturing  a vortex  in  this  sort  of  flow  by  Lin  & Rubin  D0],[11] 
and  in  some  unpublished  work  at  Douglas,  so  it  was  felt  that  the  PEBL  calcu- 
lation might  possibly  yield  a vortical  pattern  in  the  crossflow.  The  computed 
crossflow  pattern  at  x = 0.3  is  given  in  Fig.  15.  Although  the  crossflow 
does  exhibit  reversals  in  direction,  there  is  no  evidence  for  the  presence  of 
a vortex.  In  fact,  our  results  are  remarkably  similar  to  the  crossflow  pat- 
tern given  by  Wang's  Fig.  8[14],  and,  therefore,  confirm  his  conclusions  about 
the  lack  of  a vortex. 

We  have  shown  that  a comparison  of  the  3DBL  and  PEBL  calculation  pro- 
cedures produces  the  same  results  when  the  3DBL  method  is  stable.  However, 
using  the  differencing  given  in  Section  5.1,  the  3DBL  calculation  becomes 
unreliable  when  w < 0,  in  agreement  with  predictions  of  stability  calcula- 
tions. The  PEBL  procedure  is  unconditionally  stable  and  encounters  no  such 
difficulties.  Qualitative  agreement  with  the  computed  results  of  Wang  has 
been  shown  in  all  cases. 

The  deficiency  of  the  3DBL  method  used  here  can  be  overcome  by  changing 
the  differencing  to  the  zig-zag  scheme  of  Krause  [28],  which  is  orecisely 
what  Wang  used  for  his  calculations.  The  simple  box  scheme  was  used  here  to 
elucidate  its  behavior  under  adverse  conditions.  We  are  presently  investi- 
gating two  alternate  schemes,  besides  the  zig-zag  for  calculating  in  regions 
of  reverse  cross  flow.  The  inability  of  the  PEBL  procedure  to  calculate 
beyond  the  first  encountered  point  of  separation  will  need  further  study. 
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Figure  15.  Cross-stream  velocity  vectors  within  the  boundary  layer  at 


Perhaps  an  adaption  of  one  of  the  inverse  methods  used  in  two-dimensional 
boundary-layer  separation  studies  may  be  applicable  to  the  PEBL  procedure. 

Until  the  initial  separation  point,  however,  it  seems  that  the  PEBL  pro- 
cedure and  a three-dimensional  procedure  similar  to  Wang's  boundary  layer 
give  comparable  results.  Thus,  it  might  be  assumed  that  the  methods  are 
interchangeable.  However,  this  is  not  the  case.  The  zig-zag  differencing 
is  only  conditionally  stable  [28],  and  it  is  possible  that  regions  of  severe 
crossflow  might  create  instabilities  and  make  the  procedure  invalid.  Related 
behavior  is  found  when  considering  the  flow  past  conical  bodies.  Calculations 
by  similarity  variables  (to  be  shown  in  the  next  section)  cannot  progress 
further  around  the  cone  than  the  position  of  crossflow  reversal.  If  the  flow 
is  not  assumed  to  be  conically  similar,  and  is  marched  down  the  cone,  as  well 
as  around  in  the  e direction,  the  situation  is  not  altered,  i.e.,  an  inte- 
gration of  the  standard  boundary-layer  equations  fails  at  the  crossflow 
reversal  point  [29]-  Thus,  bodies  like  blunted  cones  and  elliptic  cones  can- 
not be  completely  described  except  by  using  the  PEBL  equations.  In  fact,  the 
argument  about  whether  the  lee  symmetry  plane  boundary  layer  can  be  computed 
independently  of  the  rest  of  the  flow  on  the  body,  has  not  been  settled  in  any 
general  way  for  bodies  other  than  sharp  circular  cones  [29].  For  all  cases, 
when  the  inflow  from  the  body  to  the  lee  plane  makes  the  independent  synmetry 
plane  calculation  invalid,  the  use  of  the  PEBL  equations  is  required. 

5.2  Conical  Flow 

We  have  used  the  numerical  method  of  Section  4.1  and  the  eddy-viscosity 
formulation  of  Section  2.4  to  obtain  a solution  of  the  system  of  equations 
for  conical  flow  given  by  (3.2.4)  to  (3.2.6)  and  (3.2.11)  to  (3.2.13)  subject 
to  (3.2.14).  Figure  16  shows  a comparison  of  calculated  and  experimental 
results  for  the  data  of  Rainbird  [30],  which  is  for  a 12.5-degree  half-angle 
cone  at  an  angle  of  attack  of  15.78  degrees  in  a supersonic  stream  with  free- 
stream  Mach  number  1.8.  The  Reynolds  number  of  the  cone,  based  on  its  axial 
length,  was  25  x 10^,  which  indicates  that  the  possible  effect  of  flow 
nonuniformities  caused  by  variable  transition  location  upon  the  measurements 
at  0.85L  can  be  neglected. 

The  computed  results  shown  in  Figure  16  were  made  by  using  both  the  "old" 
and  the  "new"  procedures.  According  to  the  comparison  of  results  with  experi- 
ment, we  observe  that  the  results  obtained  with  the  new  procedure  give  better 
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agreement  with  experiment  than  the  old  procedure.  For  circumferential  angles 
())  =>  0 and  45“  (see  Figs.  16a  and  16b),  however,  there  is  very  little  dif- 
ference between  the  computed  results  obtained  by  either  procedure;  as  ^ 
increases,  so  does  the  difference  as  shown  in  Figs.  16c  and  16d. 

A note  should  be  made  of  the  discrepancies  in  the  temperature  and 
velocity  profiles  at  41  = 135“.  Because  of  the  relatively  high  angle  of 
attack,  the  flow  beyond  about  41  » 120“  is  subject  to  an  adverse  pressure 
gradient,  which  eventually  leads  to  separation  defined  by  zero  shear  stress 
component  normal  to  the  generator.  Although  our  calculations  with  either 
procedure  predict  separation  at  about  4 = 161“,  close  to  the  measured  value 
of  4 = 159“,  the  measured  boundary  layer  shows  a more  rapid  thickening  than 
is  predicted  by  calculations.  This  can  be  traced  to  vortex  formation,  in 
which  case  the  ordinary  boundary-layer  equations  are  no  longer  valid  anyway. 

We  conclude  that  the  similarity  variable  proposed  by  Bradshaw  for  turbu- 
lent flows  over  tapered  wings  is  also  applicable  to  the  calculation  of  high 
Reynolds  number  turbulent  boundary  layer  on  cones  at  incidence  with  local 
similarity  assumptions. 

Although  the  improvements  obtained  by  the  "new"  procedure  may  seem  small 
in  comparison  with  the  empiricism  contained  in  eddy-viscosity  laws  and 
viscous/inviscid  interaction  effects,  they  are  nevertheless  encouraging  that 
the  local  similarity  will  work  for  other  similar  shapes  as  well.  In  principle, 
the  procedure  is  strictly  applicable  to  conical  bodies  with  arbitrary  cross- 
section  provided  that  the  radial  pressure  gradient  is  very  small. 
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APPENDIX  A 


THE  PARABOLIC-ELLIPTIC  BOUNDARY-LAYER  EQUATIONS 

The  standard  three-dimensional  boundary- layer  equations  have  a wide  range 
of  applicability.  However,  when  there  is  a region  where  the  velocity  changes 
rapidly  in  the  crossflow  direction,  a straight-forward  integration  of  the  equa- 
tions sometimes  breaks  down.  At  this  point  it  is  necessary  to  reexamine  some 
of  the  approximations  that  led  to  the  boundary-layer  equations  to  see  if  any 
have  become  invalid.  The  major  assumption  in  boundary-layer  theory  concerns 
the  relative  size  of  gradients  of  the  flow  quantities.  That  is,  gradients  in 
one  direction  are  assumed  larger  than  gradients  in  the  other  two  directions. 
However,  some  flows  do  not  fit  within  the  realm  of  boundary-layer  theory,  so 
less  stringent  assumptions  must  be  placed  on  the  Navier-Stokes  (NS)  equations. 
Flows  of  this  type  are  corner  flows  [ai].  flows  past  finite  width  flat  plates 
[A2].  rectangular  jet  flows  [23],  [A3],  flows  past  bodies  at  moderate  incidence 
[10],  [11],  [A4],  and,  as  pointed  out  in  reference  23,  wing-tip  flows.  In  these 
flows,  gradients  in  a second  direction  can  become  important,  and  must  be 
accounted  for. 

It  is  possible  to  make  a distinction  among  these  flows  which  lie  somewhere 
between  simple  boundary  layers  and  solutions  to  the  complete  NS  equations. 

In  some  cases  it  is  essential  that  the  pressure  field  be  calculated  as  an 
integral  part  of  the  solution  method.  Effectively,  this  states  that  the  scale 
lengths  in  the  gradient  directions  are  nearly  equal.  This  leads  to  the 
parabolic-elliptic  Navier-Stokes  (PENS)  equations  used  by  references  23,  A1-A4. 
However,  there  are  certain  circumstances  where  the  pressure  may  still  be  pre- 
sumed given  by  an  inviscid  flow  calculation.  This  implies  that  the  gradient 
in  the  direction  of  the  presumed  boundary  layer  is  still  the  dominant  factor, 
and  the  inclusion  of  a second  direction  only  corrects  the  flow  field  locally. 

This  leads  to  the  parabolic-elliptic  boundary-layer  (PEBL)  equations  used  by 
references  10,  11,  and  in  this  report. 

The  PEBL  equations  can  be  derived  for  the  case  of  laminar  incompressible 
flow  past  an  arbitrary  body  of  revolution  by  considering  the  Navier-Stokes 
equations  in  the  general  body-oriented  coordinate  system  discussed  in  Section 
II.  Reducing  the  general  orthogonal  equations  given  in  reference  AS  using  the 
fact  that  h,  * h^(x,()),  h^  “ h2(x,«),  h^  » 1 , yields  after  considerable  algebra 
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These  are  the  complete  Navier-Stokes  equations  for  the  body-oriented  coordi- 
nate system,  no  approximations  have  been  made. 


The  three  momentum  equations  will  now  be  nondimensional ized  in  order  to 
more  easily  determine  the  order  of  maqnitude  of  the  various  terms  in  the  equa- 
tions. The  procedure  is  the  same  as  is  used  in  standard  boundary-layer  theory 
but  with  one  important  difference.  This  appears  in  the  treatment  of  the 
cross-flow  velocity,  w.  and  the  derivatives  in  the  cross  plane,  8/36.  The 
following  nondimensional  quantities  are  defined. 


u = 


X = 


(A. 5a) 


These  are  standard  for  boundary -layer  theory  where  gradients  in  the  z-direction 

are  larger  than  those  in  the  x-direction,  i.e.,  6 <<  1 . The  subscript  r 

2 

refers  to  a suitable  reference,  and  q is  p^u  for  incompressible  flow. 

The  only  change  from  the  standard  procedure,  since  gradients  in  the  cross 
plane  are  being  considered  to  possibly  be  large,  will  be  to  flag  the  cross- 
plane gradients  with  a parameter  e as  follows: 


and  so  we  must  also  have 


(A. 5b) 


w = 


(A. 5c) 


Here  e represents  simply  a means  of  keeping  track  of  e-derivatives.  The 
value  of  this  artifice  will  appear  shortly  when  the  nondimensional  equations 
are  displayed. 

Placing  the  nondimensional  variables  (A. 5)  into  the  continuity  and 
momentum  equations,  (A.l)  - (A. 4),  rearranging  terms,  and  dropping  the  bar 

denoting  a nondimensional  quantity,  one  obtains 
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Continuity 


x-Momentum 
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Following  the  standard  procedure  of  boundary-layer  thoery,  we  now  choose 
our  y-scale  such  that 

V 1 , , 

« = rr  = R (A-10) 

r r 

Up  to  this  point,  only  the  inclusion  of  the  e flag  distinguishes  this 
development  from  that  of  ordinary  boundary-layer  theory.  If  we  were  to  set 

e = 1 and  then  drop  terms  of  0(6)  or  higher,  the  standard  three-dimensional 
boundary-layer  (3DBL)  equations  result.  However,  it  is  just  these  equations 
we  are  trying  to  modify.  They  are  to  be  modified  due  to  phenomena  in  the 
cross  plane  which  may  cause  the  e-gradients  to  become  larger  than  unit  order. 
The  terms  affected  by  these  gradients  can  be  identified  by  the  terms  flagged 
with  an  e.  The  largest  of  these  terms  in  (A. 7)  through  (A. 9)  is  of  order 

p p 

6 /e  in  each  equation.  It  is  these  extra  terms  we  wish  to  keep  in  addition 
to  the  three-dimensional  boundary-layer  terms. 

The  behavior  of  the  extra  gradients  will  not  be  uniform  throughout  the 
flow  field,  so  we  cannot  set  one  scale  for  e as  we  did  for  6.  For  example, 
near  the  windward  generator  on  a cone  no  additional  terms  are  necessary  to 
describe  the  syimietry  line  flow  other  than  the  3DBL  equations;  hence, 
e -v  0(1).  However,  near  the  lee  meridian  at  moderate  angles  of  attack, 
account  must  be  taken  of  the  inflow  to  the  syimietry  line,  and  so  a reasonable 
guess  is  that  here,  in  this  boundary  region,  e -v.  0(6).  Thus,  the  procedure 
followed  in  keeping  the  significant  terms  in  (A. 7)  through  (A. 9)  is  to 
keep  all  terms  of  0(1)  and  0(6  /e  ).  All  other  terms  are  omitted.  After 
this  has  been  done,  having  identified  the  terms  necessary  to  account  for 
cross-flow  effects,  c is  set  equal  to  one  since  its  true  scaling  is  unknown. 
The  resulting  equations  contain  terms  which  are  of  no  importance  when  e is 
truly  order  one  (hence  their  inclusion  does  not  alter  the  flow  structure), 
but  it  is  just  these  terms  which  will  make  a contribution  in  high  9-gradient 
regions.  Following  this  procedure  in  equations  (A. 6)  through  (A. 9)  yields 
the  following  set  of  parabolic-elliptic  boundary-layer  (PEBL)  equations,  where 
the  definition  of  6,  (A. 10),  has  been  used. 
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Continuity 
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All  the  terms  above  are  easily  recognized.  Those  with  no  Reynolds  num- 
ber are  the  ordinary  boundary-layer  equations,  and  remaining  terms  multiplied 
by  1/R  are  the  new  ones  needed  to  account  for  cross-plane  phenomena.  Thus, 
the  equations  are  now  elliptic  in  two  directions  (9  and  y)  but  still  allow 
a parabolic  march  in  the  x-direction.  The  pressure  is  determined  from  the 
externally  determined  inviscid  field.  This  is  by  no  means  the  same  as  having 
an  additional  equation  similar  to  (A. 13),.  instead  of  (A. 14),  had  the 
scale  length  c been  chosen  identical  to  6 at  the  outset  of  the  analysis. 
This  leads  to  the  parabolic  elliptic  Navier-Stokes  equations,  used  by  [23], 
[A3]-  [A4],  where  the  computation  of  the  pressure  is  an  essential  part  of  the 
solution.  It  is  precisely  the  calculation  of  the  pressure  that  makes  NS 
calculations  difficult.  By  choosing  the  PEBL  scaling  in  suitable  problems, 
this  difficulty  is  overcome.  This  procedure  has  been  used  previously  by  Lin 
and  Rubin  for  compressible  flow  past  nonrotating  [10]  and  rotating  [ji]  cones 
at  incidence. 
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